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Abstract 


A computational fluid dynamics (CFD) approach to predicting high-speed 
aerod 5 mamic flow fields of interest to the U.S. Army Research Laboratory (ARL) 
has been carried out. The aerodynamic problems of particular interest are 

(1) supersonic flow past the aftbody of projectiles with base mass injection, 

(2) supersonic flow past the M549 projectile, and (3) subsonic, transonic, and 
supersonic flow past an M864 projectile with base bleed and wake combustion. 

The commercially available FLUENT (Fluent, Inc. FLUENT. Version 5.1.1, 
Lebanon, NH, 1999.) CFD code was utilized. The computational effort supports 
an ongoing ARL-sponsored experimental investigation. Of particular interest in 
the present investigation is the careful characterization of the various turbulence 
models employed in the CFD code. Additionally, the ease of use and set-up as 
well as the computational time will be described. 

An experimental effort (Dutton, J. C., and A. L. Addy. 'Tluid Dynamic 
Mechanisms and Interactions Within Separated Flows." U.S. Army Research 
Office Research Grant DAAH04-93-G-0226 and the Department of Mechanical 
and Industrial Engineering, University of Illinois, Urbana-Champagne, Urbana, 
IL, August 1998.) consisting of detailed laser Doppler velocimeter (LDV), particle 
image velocimeter (PIV), and high-speed wall pressure measurements has been 
made in axisymmetric and planar subsonic and supersonic flows with embedded 
separated regions. The present work seeks to predict similar flow fields 
computationally and to address areas of agreement and disagreement. 
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1. Introduction 


A computational fluid dynamics (CFD) approach to predicting high-speed 
aerodynamic flow fields of interest to the U.S. Army Research Laboratory (ARL) 
has been carried out. The aerodynamic problems of particular interest are 
supersonic flow past the aftbody of projectiles witii base mass injection and 
supersonic flow past the M549 projectile. 

The commercially available FLUENT [1] CFD code was utilized. The 
computational effort supports an ongoing ARL-sponsored experimental 
investigation. Of partictdar interest in the present investigation is the careful 
characterization of the various turbvdence models employed in flie CFD code. 
Additionally, flie ease of use and set-up as well as the computational time will be 
described. 

Previous computational studies of base flow for transonic and supersonic free 
stream velocities carried out at ARL using Navier-Stokes techniques are reported 
in references [2-4]. Sahu and Heavey [3] compared the results of their 
computational study to experimental data and found the standard k-epsilon (k-e) 
turbulence model performed better in the wake region than an algebraic model. 

An experimental effort [5] consisting of detailed laser Doppler velocimeter 
(LDV), particle image velocimeter (PIV), and high-speed wall pressure 
measurements has been made in axisymmetric and planar subsoiuc and 
supersonic flows with embedded separated regions. The work concentrated, in 
part, on the following key issues: 

• supersonic base flow in the near wake of a cylindrical afterbody, 

• supersonic near-wake afterbody boattailing effects on axisymmetric bodies, 

• effects of rapid expansion on ihe development of compressible free shear 
layers, 

• subsonic base cavity flow-field structure, 

• base bleed with a cylindrical afterbody in supersonic flow, 

• turbtdent structures in a supersonic base flow with base bleed, 

• turbulence structure of reattaching axisymmetric free shear layers, and 

• shock-separated free shear layers. 

The present work seeks to predict similar flow fields computationally and to 
address areas of agreement and disagreement. 
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2. Supersonic Flow Past an Aftbody 


2.1 Experimental Flow Field 

The flow field investigated is a blunt cylindrical body with base bleed aligned in 
a supersonic flow (Figure 1). The supersonic free stream expands as it turns the 
comer while the turbulent boundary layer separates and then undergoes 
recompression, realignment, and redevelopment in the wake of the underbody 
[5]. Fluid from the region adjacent to the base is entrained and accelerated by the 
outer shear layer and then returned to the base region by a recompression shock 
system. This region is referred to as the primary recirculation region. 
Introducing base bleed, the primary recirculation region is moved downstream 
from the aftbody with a forward stagnation point created, dependent on the 
relative strengths of the bleed jet and the recirculating region. Experiments 
performed by several investigators (see reference [5] for a complete hst) have 
demonstrated the important effect of such a shift in the location of the primary 
recirculation region—a change in the base pressure ratio and, as a restdt, a 
change in the ctftbody drag. Base bleed then is an effective mechanism for 
reducing aftbody drag. 

The experimental flow conditions and geometry are listed in Table 1. 

2.2 Computational Grid 

Gambit [6], a single integrated preprocessor for CFD analysis, was used for 
geometry construction and import. In addition. Gambit, used for mesh 
generation, has the capability to produce structured and unstructured 
hexahedral, tetrahedral, pyramid, and prism computational cells. Mesh quality 
examination as well as boimdary zone assignment capability is also provided. 
The general sequence of operation for Gambit's use includes the following steps: 

• initial setup including solver selection, mesh size, etc., 

• geometry creation, 

• meshing including edge and boundary layer meshing followed by face 
and/or volume meshing, 

• mesh examination, 

• zone assignment, and 

• mesh export 
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For the present investigation, the flow field is modeled as axisymmetric flow past 
a cylindrical aftbody wifli no swirl. Thus, title generation of suitable grid 
geometry was straightforward. An appropriate length of the cylinder embedded 
in the flow field was calculated using botmdary layer flieory corresponding to 
existing experimental data. For all grids generated in the present work, a length 
of cylinder equal to 0.14 m was chosen. 

An important aspect of the present work was to determine the effect of mesh 
density on the calculated results. For the grid face size area of 0.4 m x 0.6 m, the 
number of cells varied from a low of approximately 7000 to a maximum of nearly 
64,000, thus affording an order of magnitude difference in this potentially 
important parameter. 

The actual grids used in this investigation are shown in Figure 2. Modifications 
to the simple cylindrical grid were incorporated in order to model a cylindrical 
boattail mounted behind the cylinder. The boattail grid geometry is presented in 
Figure 3 and is based on the experimental prototype of Herrin and Dutton [7]. 
The boattail chosen for this experimental investigation has a conical shape with 
an angle relative to the horizontal of 5° and length 31.75 mm (0.5 cal.). 
Reid and Hastings [8] have shown that the optimal boattail shape is essentially 
conical at moderate supersonic speeds for typical boattail lengths. The 5° angle 
has also been shown to be near tihe optimal angle from previous investigations 
[ 8 ]. 


2.3 Computational Models 

Three different turbulence models [9] are used in the present investigation: 
(1) the Spalart-AUmaras one-equation model, (2) the standard k-£ two-equation 
model, and (3) the Re 3 molds stress five-equation model (Table 2). 

In turbulence models that employ the Boussinesq approach, the central issue is 
how the eddy viscosity is computed. The model proposed by Spalart and 
AUmaras solves a transport equation for a quantity that is a modified form of the 
turbulent kinematic viscosity. 

The standard k-e model is a semi-empirical model based on model transport 
equations for the turbulent kinetic energy (k) and its dissipation rate (e). The 
model transport equation for k is derived from the exact equation, while the 
model transport equation for e was obtained using physical reasoning and bears 
little resemblance to its mathematically exact coxmterpart. 

In the derivation of the k-£ model, it was assumed that the flow is fully turbulent, 
and the effects of molecular viscosity are negligible. The standard k-£ model is 
tiierefore valid only for fuUy turbulent flows. 

The Reynolds stress model (RSM) is the most elaborate turbulence model that 
FLUENT [1] provides. Abandoning the isotropic eddy-viscosity hypottvesis, the 
RSM closes the Reynolds-averaged Navier-Stokes equations by solving transport 
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equations for the Reynolds stresses, together with an equation for the dissipation 
rate. This means that four additional transport equations are required in two- 
dimensional (2-D) flows, and seven additional transport equations must be 
solved in three-dimension (3-D). 

2.4 Comparison With Varying Grid Densities 

Results are described for various flow field properties with the turbulence 
closure model kept constant but with the grid density or number of cells within 
the grid volume varying. The turbulence models and the corresponding number 
of cells are listed in Table 3. 

In this report, niimerical results for the RSM are presented. Streamwise mean 
velocity contours are presented in Figure 4 and the turbulent kinetic energy 
contours are shown in Figure 5. These results are for a mass bleed ratio equal to 
0.01. For the case of the streamwise velocity, the contours are nearly identical for 
the approximate order-of-magnitude difference in grid density. The results are 
quite different for the turbulent kinetic energy contours, particiilarly near the 
base of the aftbody. The region of relatively low turbulent kinetic energy extends 
much fcirther downstream from the base region indicating a much diminished 
momentum mixing or transport occurring from the outer free stream towards the 
wake. 

Though not included here, results for the streamwise mean velocity contours for 
both the Spalart-AUmaras and standard k-e models behaved similarly. That is, 
the contours remained nearly identical from coarsest to finest grid density. 

2.5 Comparison With Varying Ttirbulence Models 

In order to compare the results for the varying turbulence closure models, 
contour plots of streamwise mean velocities are presented. The profiles are 
shown for the finest grid density. 

In Figure 6, the steamwise velocity contours are compared for two closure 
models. Note that a similar observation can be made about the region near the 
base as was made for the turbxilent field for the case of varying grid density. The 
more detailed and presumably higher order model, the RSM, tends to show an 
elongated wake region near the base as compared to the Spalart-AUmaras model. 
Thus, increasing the complexity of the model seems to result in the same effect as 
increasing the grid ceU density. Though not shown here, the radial velocity 
gradient is greatest for the case of the Spalart-AUmaras model (i.e., the lowest 
order model) and least for the RSM (i.e., highest order model). SimUarly, static 
pressure contours exhibit higher gradients for the Spalart-AUmaras models as 
compared to the RSM. 
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2.6 Comparison With Experimental Data 

Comparisons are made with experimental data for supersonic flow past a 
cylindrical aftbody without base bleed, with base bleed, and for a cylindrical 
aftbody plus boattail with and without base bleed [5]. For all the comparisons 
shown, the Reynolds stress turbulence closure model is used as well as the finest 
ceU density. For the cases involving base bleed, the flow exits the cylindrical 
aftbody parallel to the free stream direction with the base bleed mass injection 
rate equal to 0.01 times the mass flux injection rate if the entire base area served 
as the exit area for the nozzle. The velocity of the base bleed is considered to be 
constant over the exit plane. 

2.6.1 No Base Bleed (I = 0.00) 

In Figure 7, streamwise mean velocity profiles are presented for the case of zero 
base bleed [5]. The non-dimensionalLzed streamwise velocity, U/Ufs, is plotted 
vs. non-dimensionalized radial location, r/R, for three downstream locations, 
x/D = 1.26, 1.42, and 1.73. The numerical results are compared to flie 
experimental data. The agreement is quite good but with the numerical model 
consistently overestimating the extent of the mean wake region and thus 
overestimating the magnitude of the mean, axial velocity gradient. 

In Figure 8, non-dimensionalized turbulent shear stress (-uv/Ufs‘^2) profiles are 
presented at the same downstream locations (x/D = 1.26,1.42, and 1.73). The 
munerical results imderestimate the maximum value of the turbulent stress and 
underestimate the extent of the turbulent velocity field as compared to 
experimental data. This is in direct contrast to the results for the mean velocity 
field. 

In Figure 9, the base pressure (here nondimensionalized by free stream static 
presstue) distribution is compared to experimental results. Averaging over the 
length of the base, the pressure ratios are nearly identical. 

2.6.2 Base Bleed (I = 0.01) 

In Figure 10, streamwise mean velocity profiles are shown for foxu: downstream 
locations (x/D = 0.95,1.26,1.95, 2.04). The results are much closer to ihe actual 
experimental data [8] suggesting that the numerical model more accurately 
predicts the important flow features in Ihe case of base bleed. 

In Figure 11, radial mean velocity profiles are presented for x/D = 0.95,1.26,1.95, 
and 2.04. The agreement with experimental data is not as good as is the case for 
the streamwise velocity data. Consistently, the numerical data overestimates the 
magnitude of the radial mean velocity at each downstream location. This would 
indicate there is more turning in the flow field numerically than measiued 
experimentally. 
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In Figure 12, the turbulent kinetic energy profiles (1/2 [u^ + + w2]i/2) at the 

downstream locations x/D = 0.95, 1.26, 1.95, and 2.04 are presented. The 
agreement both in magnitude of the turbulent kinetic energy and in location is 
quite good. 

In Figure 13, the turbulent stress, -uv/Ufs''2, profiles at x/D = 0.95, 1.26, 1.95, 
and 2.04 also demonstrate excellent agreement with the experimental data. 

In Figure 14, the downstream development of the axial mean velocity is 
presented. The magnitudes are in quite close agreement with the experimental 
data as is the location of the zero velocity or leading edge of the recirculation 
region. 

In Figure 15, file base pressure, Pb/Pfs, is shown for the region outside the base 
bleed exit with the agreement with experimental data quite good. 

2.6.3 Aftbody Boattail Without and With Base Bleed (I = 0.00 and I = 0.01) 

Axial mean velocity, U/Ufs, and turbulent kinetic energy, k/Ufs'^2, as a function 
of downstream location are shown in Figures 16 and 17, respectively, for the no 
bleed and bleed cases [5]. Agreement is quite good throughout the range of 
downstream locations though does improve with distance from the exit plane. 
This is especially true downstream of the leading edge of the recirculation 
region. 


3. Supersonic Flow Past M549 Projectile 


A series of computations have been performed for the M549 configuration shown 
in Figure 18. This portion of the present investigation was performed in order to 
gain experience using the grid generation technique for a more complicated 
geometry than the previously investigated flow past an aftbody. The Mach 
number was set equal to 2.47 for all computations at zero angle of attack. The 
computational flow conditions and geometry are listed in Table 4. 

Qualitative features of the computed flow fields are presented in the form of 
contours of the axial and radial mean velocity fields, and contours of the 
turbulent kinetic energy and uv-Reynolds stress. 

The grid geometry for the flow past the M549 projectile (not shown) was 
obtained using the grid-adaption feature of the FLUENT [1] code to obtain a 
higher density of cells in region of high static pressure gradient. No 
experimental data are available for comparison so only numerical results are 
provided. 

Contours for the axial velocity and radial mean velocities are presented in 
Figures 18 and 19. In Figure 20, the contour plot for the mean static pressure 
field is shown. The turbulent velocity field is detailed with contour plots of the 
turbulent kinetic energy and the Reynolds shear stress in Figures 21 and 22, 
respectively. 
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The axial mean velocity contours shown in Figure 18 demonstrate the important 
characteristics of the flow field. The flow first encoimters a compression region 
near the nose of the projectile, followed by a weaker compression at the waist, 
and then subsequently followed by a turbulent wake region at the base. 
Evidence of the turbulent botmdary layer growth along the length of the 
projectile is evident. The radial mean velocity contours presented in Figure 19 
indicate the existence of strong recirculation behind the base region and large 
influx of the outer free stream flow. The mean static pressxue contotirs 
(Figure 20) reinforce the earlier observations concerning tine existence of 
compression regions near the nose, near the waist, and in the wake region behind 
the base. 

The turbulent flow field contotus presented in Figures 21 and 22 add further to 
an understanding of the resultant flow. The uv-component of the Reynolds 
stress is selected as it provides the greatest insight into the origin of the 
turbulence. Note that die regions of highest turbulent kinetic energy correspond 
to regions of relatively low magnitude of the shear stress term. Conversely, the 
maximum values for the shear stress occur immediately behind the base where 
die turbulent kinetic energy remains at a low level. Thus, the turbulent flow field 
begins to develop or is generated close to the base and then, as the flow convects 
downstream, reaches a maximum before being attenuated by the effect of the 
viscosity of the fluid. 


4. Subsonic and Supersonic Flow Past M864 Projectile 


The M864 is an Army shell that bums a solid propellant diat exhausts into the 
base region to achieve reduced aerodynamic drag. It is a 155-mm, extended 
range, cargo-carrying projectile. A schematic of the M864 is shown in Figure 23. 
Note that the boattaded section includes a domed cavity and a combustion 
chamber for the solid propellant. The products of the combustion exit through 
an orifice plate rather than a traditional DeLaval (C-D) nozzle, as the M864 is 
considered to be of base bleed design rather than a rocket-assisted projectile. 

An initial numerical investigation by Nietubicz and Gibeling [10] modeled the 
solid propellant combustion using the approach developed by Giebling and 
Buggeln [11]. This study [10] modeled the combustion as injectant gas mixtures 
of H 2 and H 2 /CO for a base bleed ratio of I = 0.0022. These combustion results 
were compared to results for hot air bleed and zero bleed conditions. The 
computations were carried out for M = 2.0 and zero angle of attack. The results 
indicated that the effect of the hot mass injection was confined to a region very 
near the base, whereas the combustion products affected a much wider region of 
the base region flow field. 
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4.1 Wake Combustion Model 


The present numerical investigation of the M864 includes the effects of the 
inclusion of the domed cavity, the injection of high temperature fluid at the base 
exit plane and the effect of wake combustion. The wake combustion has been 
incorporated into the calculation based upon the published FLUENT [1] tutorial 
for nozzle flow for a solid-propellant rocket. A brief description of this model is 
now provided. 

The mass flow of the solid propellant exiting a nozzle, or, in the present case, an 
orifice, is given by: 

l = 5r/dt* Apb = Ti ♦ [Pc] TN, (1) 

where, 

dr/dt = surface burning rate, 

A = base area, 

pb = density of solid propellant, 

Ti = empirically determined constant. 

Pc = fluid pressure in the combustion chamber, and 
N = empirically determined constant. 

For the present work, the mass flow rate was set at 0.01, and then the empirical 
constants were calculated at the appropriate flow conditions to obtain this mass 
flow rate. 

Computations have been carried out for the case of zero angle of attack for a free 
stream Mach number range of 0.8-3.0 (Table 5). Qualitative features of the 
computed flow fields are presented in the contours of the axial and radial mean 
velocity fields, contours of the turbulent kinetic energy and uv-Reynolds stress, 
and contours of the mean static temperature fields. C^antitative descriptions of 
the various flow fields are also provided including the downstream or 
streamwise development of the mean velocity distributions. 

For the present investigation, the RSM for turbulence was used exclusively. This 
particular turbidence model was chosen based on the earlier work of the 
principal investigator focusing on the simpler aftbody geometry, and the M549 
projectile. 

4.2 Discussion of Results 

Contours for the axial mean velocity fields for the three Mach numbers with base 
bleed and both with and without wake combustion are presented in 
Figures 24-26. For a given mass bleed, increasing the free stream Mach number 
increases the turning of the flow around tire base exit plane. This is similar to the 
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effect of increasing the relative velocity differences in mixing layers and m 
coaxial jets. Greater mean velocity ratios serve to increase the mixing more 
rapidly. 

Contours for the radial mean velocity fields for the supersonic Mach numbers 
with base bleed and both with and without wake combustion are presented in 
Figures 27 and 28. Increasing the velocity ratio puUs tiie standing vortex located 
near the base region toward the rear of the projectile. The wake combustion in 
turn flattejis the vortex and tends to diffuse or flatten the area of steepest radial 
velocity gradient toward the outer reaches of the base region. 

Contours for the turbulent kinetic energy fields for the three Mach numbers with 
base bleed and both widi and without wake combustion are presented in Figures 
29-31. The presence of wake combustion seems to aid in the transfer of the flow 
energy from the mean velocity field to the turbulent velocity field. That is, the 
regions of highest turbident kinetic energy have convected downstream more 
rapidly with the inclusion of the wake combustion. However, it is worth noting 
that the high supersonic case seems to be least affected. 

Contours for tiie turbulent uv-Reynolds shear stress fields for the three Mach 
nximbers with base bleed and both with and without wake combustion are 
presented in Figures 32-34. Increasing the free stream Mach ntimber tends to 
increase the magnitude of the uv component of the Reynolds shear stress and 
pull that region of maximum value closer to the base. Similarly, the introduction 
of wake combustion also increases the magnitude of ttie uv component of 
Reynolds stress and pull regions of highest gradient toward the base. 

In Figures 35-38, contour plots for M=3.0 and with wake combustion are 
presented. The contotu plots indicate the important features of the resultant 
flow field and are not provided as a smnmary of the computational results. 


5. Summary 


Computations have been performed for flow past an aftbody, an aftbody with 
boattaUing, and the M549 projectile configuration. 

A numerical investigation has been made for supersonic flow past a cylindrical 
aftbody with and without base bleed and with and without a cylindrical boattail. 
The effects of grid cell density and turbulence model were each examined. The 
computational fluid dynamics software FLUENT [1] Version 5.1.1 was used for 
all calculations. The presentation of the graphs extensively used Adobe 
Illustrator Version 7.0. The mean and turbulent velocity and pressure fields all 
were documented. 
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The grid density was found to have its most significant effect on the calculations 
of the turbulent velocity field while the mean velocity field was essentially 
independent of grid cell density. 

Three different turbulent models were employed in the present investigation, the 
Spalart-AIlmaras one-equation model, die standard k-e two-equation, model and 
the Reynolds stress five-equation model. The more detailed and presumably 
higher order model, the RSM, tends to show an elongated wake region near the 
base as compared to the other two models. Increasing the complexity of the 
model seems to result in the same effect as increasing the grid cell density. 

The numerical calculations were found to be closest to die experimental data for 
the case of mass base bleed than in the case of no mass base bleed. The presence 
or absence of the aft boattail did not seem to have an effect on the comparison of 
numerical and experimental data. 

For the flow past the M549 projectile, the mean and turbulent velocity fields were 
described with experience gained in modeling flows past finite bodies of 
revolution. Before this portion of the investigation, only sections of the projectile, 
the aftbody, were investigated numerically. 

Computations for the reacting M864 base bleed projectile provided results for the 
effect of temperature on the base region standoff distance that were not fuUy 
consistent with previous results of Nietubicz [10], which modeled the 
combustion process as a finite-rate reacting gas mixture. However, there are 
significant differences in the two approaches with respect to the combustion 
model, base bleed mass rate, numerical scheme, and tiubulence model. 
Additional effort is required to identify the source of the discrepancies and 
establish the optimum modeling technique for this complex flow field. 
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Figure 1. Schematic diagram of supersonic base flow. 


11 














































Figure 2. Grid geometiy and variation of grid densities for flow past aftbody. 




[a] [b] 

Figure 4. Comparison of streamwise mean velocity contours for high and low grid 
number densities ([a] high grid number density = 64k vs. [b] low grid number 
density = 9k). 












Figure 6. Comparison of streamwise mean velocity contours for varying turbulence 
models. 
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Figure 7. Streamwise mean velocity (U/Ufs) profiles as a function of radial position 
(r/R) for various downstream locations (compared to experimental data) for no 
base bleed (I = 0.0). 
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Figure 8. Turbulent shear stress (-uv/Ufs^2) as a function of radial position (r/R) for 
various downstream locations (compared to experimental data) for no base 
bleed (I = 0.0). 






r/R, Radial Location 



Pb /Pfs, Base Pressure Distribution 

Figure 9. Non-dimensionalized base pressure distribution, Pb/Pfs, vs. radial location on 
base for no base bleed (I = 0.0) with comparison to experimental data. 


17 



r / R, Radial Location r / R, Radial Location 



Figure 10. Streamwise mean velocity (U/Ufs) profiles as a function of radial position 
(r/R) for various downstream locations (compared to experimental data) for 
base bleed (I = 0.01). 
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Figure 12. Turbulent kinetic energy (k/Ufs^2) as a function of radial position (r/R) for 
various downstream locatioi\s (compared to experimental data) for base bleed 
(I = 0.01). 


20 









r/R, Radial Location 













0 


2 


4 6 8 

X / R, Downstream Location 


Figure 14. Streamwise mean velocity, U/Ufs, vs. downstream location x/R, for base 
bleed (I = 0.01) with comparison to experimental data. 



Figure 15. Base pressure profile for base bleed (I = 0.01). 
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x/ R, Downstream Location 


Figxire 16. Streamwise mean velocity, U/Ufs, as a function of downstream location, x/R, 
for no base bleed and base bleed (I = 0.01) with comparison to experimental 
data. 



x/R, Downstream Location 


Figure 17. Turbulent kinetic energy, k/Ufs^2, vs. downstream location, x/R, for no bleed 
and base bleed (1 = 0.01) with comparison to experimental data. 
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Figure 18. Contours of mean axial velocity for M549 projectile. 
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Figure 20. Contours of static pressure for M549 projectile. 
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Figure 21. Contours of turbulent kinetic energy for M549 projectile. 



Figure 22. Contours of turbulent shear stress for M549 projectile. 






Figure 23. Schematic of M864 projectile. 
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Contours of mean axial velocity for M = 0.8 without and with wake 
combustion. 
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combustion. 
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Figure 35. Contours of Mach number for M = 3.0 with wake combustion. 
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Figure 36. Contours of axial mean velocity for M = 3.0 with wake combustion. 














Table 1. Experimental flow conditions and geometry. 


Flow Property and Geometry 

Magnitude 

Free Stream Static Pressure 

28700 Pa 

Approach Free Stream Mach Number 

2.47 

Tunnel Stagnation Temperature 

300 K 

Bleed Flow Mass Flow Rate Ratio 

0.01 

Base Radius 

31.75 mm 

Bleed Orifice Radius 

12.7 mm 

Bleed Flow Stagnation Temperature 

300 K 

Tunnel Stagnation Pressure 

470000 Pa 


Table 2. Turbulence closure models. 


Tubulence Closure Models 

Key Features/Assumptions 

Spalart-Allmaras Model 

Boussinesq approximation, 1-equation 

Standard k-c Model 

Relates dissipation rate and turbulent 
kinetic energy, 2-equation 

RSM 

Relates Reynolds stresses and dissipation 
rate, 5-equation 


Table 3. Turbulence closure models and cell density. 


Turbulence Model 

Coarse Density 

Medium Density 

Fine Density 

Spalart-Allmaras 

7197 

15724 

63809 

k-e 

8770 

15724 

63809 

RSM 

8690 

15724 

63809 
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Table 4. Computational flow conditions for M549 projectile. 


Flow Property and Geometry 

Magnitude/ Characteristic 

Projectile 

M549 

Free Stream Static Pressure 

28700 Pa 

Approach Free Stream Mach Number 

2.47 

Tunnel Stagnation Temperature 

300 K 

Turbulence Model 

RSM 5-equation 

Cell Density 

71234 

Mass Base Bleed 

NA 

Tunnel Stagnation Pressure 

470000 Pa 


Table 5. Computational flow conditions for M864 projectile. 


Flow Property and Geometry 

00 

d 

11 

M = 1.5 

O 

cd 

II 

Projectile 

M864 

M864 

M864 

Free-Stream Static Pressure (Pascals) 

28700 

28700 

28700 

Stagnation Temperature (Kelvin) 

300 

300 

300 

Cell Density 

105300 

105,500 

105300 

Stagnation Pressure (Pascals) 

43738.8 

105357.7 

1054237 

Turbulence Model (RSM/5-equation) 

RSM 

5-eq 

RSM 

5-eq 

RSM 

5-eq 

Mass Base Bleed (I) 

0.01 

0.01 

0.01 

0.01 

0.01 

0.01 

Wake Combustion (No/Yes) 

N 

Y 

N 

Y 

N 

B 
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